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BACKGROUND: In this study, the angle of the myocardium's trajectory in three dimensions (O) was estimated by simultaneous 
use of long-axis and short-axis views of left ventricle septum two-dimensional images. Then correlation of three-dimensional 
trajectory's angle with the rotation angle from the long (X) and short (0) axis views was estimated and compared at the three levels 
of base, mid and apex of the interventricular septum wall. 

METHODS: Two-dimensional echocardiography images of long- and short-axis views of 19 healthy men were recorded and 
analyzed. Using an electrocardiogram of each individual, the images of the two views were synchronized. The interventricular 
septum wall motion at the three levels of base, mid and apex were estimated, using a block matching algorithm throughout 
three cardiac cycles. Considering the defined system of coordinates and the position vectors in long and short-axis views, the 
3 -dimensional angle of the trajectory was calculated. 

RESULTS: Maxima of the O, 0, and X angles were extracted at 16.33 + 3.01, 10.61 ± 3.38, and 15.11 ± 3.30 degrees at base level, 
22.77 ± 4.95, 7.78 ± 2.96, and 16.72 ± 2.66 degrees at mid level and 14.60 ± 5.81, 10.37 ± 5.48, and 8.79 ± 3.32 degrees at apex 
level, respectively, of the septum wall, respectively. This study shows significant correlation between the angle of 3-dimensional 
trajectory (O) with the angle in short axis view (0) of the septum wall at the apex level; and also with the angle in long axis view 
(X) of the septum wall at base and mid levels. 

CONCLUSION: Due to the motion of the wall of the left ventricle in three dimensions, and the non-isotropic structure of 
myofibers, the angle of 3-dimensional trajectory was estimated using the speckle tracking method of 2-dimentional 
echocardiography images. 
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Introduction 

General performance parameters of the left ventricle (LV), 
such as the ejection fraction, cannot differentiate the local 
changes caused by different physiological conditions. Local 
study of the left ventricular performance was carried out by 
quantification of the myocardium trajectory, using the param- 
eters of strain, strain rate and stress. The helix structure of 
the myofibers of myocardium and the change of these myofi- 
bers' angle along the wall thickness, lead to torsional motion 
of the ventricle (apex rotational motion in the opposite direc- 



tion of the base).^^^^ According to the Torrent-Gausp model,^^ 
the lefi: ventricular tissue consists of two descending (clockwise) 
and ascending (counter-clockwise) branches which form the 
base section at first, then descend spirally, and ascend to the aor- 
ta. In addition to the uniform distribution of stress and shorten- 
ing along the fibers, torsion also plays a major role in the ejec- 
tion fraction of the heart. So with a shortening of myofibers of 
about only 15—20%, the heart's ejection fraction is estimated at 
50—60%.^^^^ Any torsional disorders increase the stress and 
strain in the inner layer, hence increasing the demand for oxy- 
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gen, and ultimately decreasing the systolic performance of the 
LY7)8) 'p]^^3^ potential characteristic of measuring LV torsion 
may bean early diagnosis of myocardial damage.^^^^ 

The quantification of such motion has been done using the 
angle of rotation, torsion angle and circumferential-longitudi- 
nal shear angle^°^ with invasive implantation of tantalum mark- 
ers and X-ray imaging,^^^ and non-invasive techniques such as 
magnetic resonance imaging (MRI) and echocardiography ^^^^^^ 
The measurement of torsion in echocardiography images was 
initially done using the tissue Doppler technique. The latest 
technique for the measurement of the torsion angle is using 
the block matching algorithm to track sequential echocar- 
diography images,^^^^^^ which does not have the problems as- 
sociated with the tissue Doppler technique. Despite the 3 -di- 
mensional motion of the myocardium, this angle has so far 
been estimated in echocardiography images by eliminating 
the longitudinal motion in 2 -dimensional short-axis view im- 
ages,^ ^^^^^^^^ and by eliminating the circumferential motion in 

2- dimensional longitudinal view images. Despite many stud- 
ies about the left ventricular angle of rotation on a two-dimen- 
sional scale, there are as yet, no studies about the introduction 
and extraction of the torsion angle of myocardial levels on a 

3 - dimensional scale. 

Thus in the current study, 2-dimensional echocardiography 
images are first synchronized for two short-axis and longitu- 
dinal-axis views, using the view's electrocardiogram (ECG) 
signals. The sequential echocardiography images are then 
processed, using the block matching algorithm. Coordina- 
tion of the interventricular septum trajectory will be estimat- 
ed at the three levels of base, mid and apex. Since the quality 
of images is very important for the accuracy of tracking, the 
highest quality images of the base, mid and apex levels of the 
septal wall are used. After extracting the 3 -dimensional tra- 
jectory, the angle of the left ventricular 3 -dimensional trajec- 
tory will be estimated, considering 3 -dimensional trajectory 
vectors. The correlation between the 3 -dimensional angle of 
trajectory in the base, mid and apex levels of the interventric- 
ularseptal wall, and extracted rotation angles from both the 
short-axis and longitudinal axis views, will be analyzed and 
assessed. 

METHODS 

Subjects 

Nineteen adult males (46 ± 6 years old), who based on 
echocardiography examinations, had a healthy heart, volun- 
teered to participate in this study. Based on Framingham 
study indices,^^^ including a person's gender, age, blood pres- 
sure, smoking, blood cholesterol levels, having diabetes and 
hypertrophy, the risk of developing coronary artery disease in 
the volunteered participants was less than 10%. Each healthy 
individual underwent normal physical examinations, includ- 
ing blood pressure measurements, echocardiographic general 



parameters studies, and recording of ECG signals, at the same 
time as echocardiography imaging. The volunteered individu- 
als had no history of cardiovascular diseases, angina, high 
blood pressure, or diabetes, and were not taking any medica- 
tion. Before echocardiography imaging, the blood pressure of 
the individuals was measured in the supine position, using a 
Microfile_BP 3BU (Microfile, Heerbrugg, Switzerland) cali- 
brated oscillometric sphygmomanometer with a blocking cuff 
(± 1 mmHg, 30—280 mmHg). General performance parame- 
ters of the heart, including linear measurements of ventricular 
dimensions and the left atrium of heart, were then extracted 
in longitudinal parasternal and apical four-chambers views, 
using M-mode and B-mode ultrasound imaging techniques. 
The ventricular diastolic performance was also extracted, us- 
ing blood flow parameters of the mitral valve, with spectral 
Doppler imaging technique. Before the test, informed partici- 
pation consent was received from all the volunteers. This pro- 
tocol was approved by the medical ethics committee of Tarbi- 
at Modares University. 

Imaging and extracting the myocardial muscle 
movement stage 

After a general examination of the subjects, tissue harmonic 
imaging (transmitted frequency of 1.7 MHz and received fre- 
quency of 3.4 MHz) was done using an echocardiography sys- 
tem (Philips iE33 xMATRlX, Best, Netherlands) with elec- 
tronic ultrasonic transducer (S5-1), and 75-95 frames per 
second in Digital Imaging and Communications in Medicine 
format, a spatial resolution of 0.26 ± 0.01 mm, and a maxi- 
mum depth of study of 12 to 16 cm. Echocardiography imag- 
es were recorded from the longitudinal four chambers and 
short-axis views at the three apex, mid and base levels accord- 
ing to the standards of the American Heart Association 
(AHA).^^^ To remove artifacts caused by movement of the dia- 
phragm, the volunteers were asked to hold their breath during 
the recording of three cardiac cycles of echocardiography im- 
ages. After transferring the sequential echocardiography im- 
ages to the computer, with a minimum temporal resolution of 
14 ms, the mid-ventricular septal wall was tracked based on 
the AHA sixteen segments standard model,^^^ using the block 
matching algorithm, in three cardiac cycles at the three base, 
mid and apex levels, and the average of these three measure- 
ments, with an error of less than 10%, was considered as the 
output. Low quality images that had artifacts, due to the 
movement of diaphragm or the subjects during the imaging 
process, were removed from the study. 

In the block matching algorithm, after reading images of a 
cardiac cycle as the input, the region of interest (ROI) in the 
source images was then selected by the user. Considering the 
size of the ROI and its location in the source image, the search 
area was defined in dimensions twice as large as the area of the 
ROI, centered on the location of the ROI in the target image. 
In the next stage, the sum of absolute difference (SAD) quan- 
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tity between the ROI, and all blocks with the same size was 
calculated in the search area. The location of the block with 
the minimum SAD was selected as the location of the ROI in 
the target image/^^^^^ 

According to the flowchart shown in Fig. 1, the ROI is se- 
lected by the user in the first frame, and the block correspond- 
ing to the minimum SAD is found in the second frame (Fig. 2). 
In order to select the corresponding block in the third frame, a 
block that has the minimum SAD, compared to the block that 
was found in the second frame, would be selected and so forth. 

After the execution of the block matching algorithm, con- 
sidering the geometric model of the heart as an truncated el- 
lipsoid (Fig. 3A),^^ and the geometry of the taken views, the 

coordinates (Ximage, Yimage) (Fig. 3B) and (_LZimage, Zimage) (Fig. 

3C) were achieved, relative to the image coordinate system. 
The reference coordinate system was considered as shown in 
Fig. 3 by selecting the center of the mid-level of the LV as the 
origin of the coordinates and choosing the orientation of its 
axes. The center of ventricle (Xo, Yo) (Fig. 3B) was identified 
in short-axis view on coordinates of the image using Fiji soft- 
ware (http://fiji.sc) (the image processing software package for 
bioscience researchers), by fit of circular or oval shapes on seg- 
ments of the LV anatomical images. The origin of coordinates 
(_LZo, Zo) (Fig. 3C) was identified in longitudinal axis view, 
using the anatomical papillary muscle marker. At the end of 
this stage, the coordinates of the myocardium's septum were 
transferred, from the coordinates system defined in the image, 
to the reference coordinate system (equations 1^). 



Fig. 1. Process view of the block matching program. ROI: region of 
interest, SAD: sum of absolute difference. 
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Fig. 2. Implementation of block matching algorithm. Region of interest (ROI) is selected by the user in the first image and the 
second image corresponding block is the minimum sum of absolute difference. Sample ROI search (frame 78) (A) and destination 
(frame 112) (B) in the septal segment echocardiographic images based on the 4-chamber view is specified. 
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Fig. 3. A: Truncated ellipsoidal model and the reference coordinate system. B: Standard short-axis view and the reference coordinate system [Mid 
Septum (Sept), Mid Inferior (Inf), Mid Posterior (Post), Mid Lateral (Lat), Mid Anterior (Ant), Mid Anterior Septum (Ant Sept)]. C: Standard long axis view 
and the reference coordinate system [Apical Septum (Sept), Mid Septum (Sept), Basal Septum (Sept), Apical Lateral (Lat), Mid Lateral (Lat), Basal 
Lateral (Lat)]. 



Synchronization of the LV views in both 
longitudinal and lateral axes 

Longitudinal and short-axis views were imaged with a 2—3 
min delay from the heart of each individual. Due to the 
changes of heart rate during successive cycles of the heart and 
unequal frame rates in both views in this study, combining 
the images in short-axis view with the ones taken in longitu- 
dinal axis view needed examinination of the synchronization 
of the short- and long-axis views data using the ECG signal. 
Therefore, the ECG signal was extracted from the echocar- 
diography images using this signal in such images and color 
images processing. In order to extract this signal, the change 
in coordinates of the last colored pixel of the ECG signal in 
the images of a cardiac cycle was used. To obtain the coordi- 
nates of this pixel, all the images of a cycle were first limited 
to the range of the ECG signal and after separating the color 
channels of the images, the image of the channel with highest 
intensity was selected (here G). Coordinates of the last pixel 
were recorded using the number 1 in all the images, after the 
conversion of the main channel images to binary data. The 
ECG signal was reconstructed by transferring these values to 
the Excel software. Signals recorded in long and short-axis 
views were compared visually afiier the extraction of the ECG 
signal in both these views of echocardiography images. The 
differences were observed in at the end of the diastolic phase 
for most of the images. In images where such differences were 
high, the repeated frames in different phases of the cardiac cy- 
cle were removed due to inequality in the size of such differ- 
ences for different phases of ECG signal in order to synchro- 
nize the data. The output data in this stage were used to 
calculate the short-axis, long-axis angles and angle of 3 -di- 
mensional trajectory. 

CALCULATING ANGLES 

Considering the defined system of coordinates, a = (a^z, az) 
and J = (^iz, bz) are position vectors in long-axis view. Thus, 
the rotation angle in long-axis view D((t)} is obtained from 



equation 5 : 
5) = s'm ' (^^,) = sin 



|2| • % ~ ^sqrt {ai + a^z) • sqrt {hz + h^z) ^ 



a = (ax, ay) and I = by) are position vectors in short-axis 
view. Thus, the rotation angle in short-axis view {6(t)} is ob- 
tained from equation 6: 

\a X l\ . _^ r sqn {{ax - by- ay - bxf] ^ 



6) m = s'm' (^^,) = sin 



sqrt {a} + a^) • sqrt {hx + h^) 



a = (ax, ay, az) and I = by, i?z) are position vectors in short- 
axis view. Thus, the 3-dimensional angle of trajectory {0(t)} is 
obtained from equation 7: 



'hz-az' byf + {ax -hz-az- + -by -ay &)^} . 
sqrt {ax + a^ + a?) • sqrt (bx + + hi) ^ 



The average of angle changes for three cardiac cycles was 
then extracted as the output of this stage. In order to perform 
statistical comparison of the results, the maximum angle of 
3-dimensional trajectory (Om), and maximum angles of rota- 
tion in short-axis (0m) and long-axis (Xm) were reported as the 
average ± standard deviation. Pearson's analysis of correlation 
with a significance level of less than 0.05 was used to study 
the correlation between 0m, Xm and Om. The sample size was 
estimated with a confidence level of 95% and power of test of 
90% in each group. The statistical analysis was done using 
SPSS 16 software (SPSS Inc., Chicago, IL, USA). 

RESULTS 

The demographic and echocardiographic measures of 19 
men with an average age of 46 ± 5 years are given Table 1. 
The general parameters extracted from Table 1 shows the 
healthiness of heart in individuals of this study. 

Recorded ECG signals in the short and long axis views im- 
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ages were visually compared for each individual, during three 
cardiac cycles. Due to the non-synchronization of three cardiac 
cycles in two recorded views, repeated data in different phases 
of the cardiac cycle were removed. Fig. 4 A was drawn before 
synchronization, and Fig. 4B was drawn after synchronization 
of the images, in both short and longitudinal axes views. 

Fig. 5 shows the changes in location of the septum segment 
in base (Fig. 5A), mid (Fig. 5B), and apex (Fig. 5C), separately, 
for a sample individual, obtained by averaging the data of three 
cardiac cycles. The location of the septal wall is presented in 
three base, mid and apex levels, using x, y, and z coordinates rel- 
ative to the origin, for the mid-level of the LV in pixels. 

The changes of x, y coordinates in three base, mid and apex 
levels of the septal wall are consistent with ventricular contrac- 
tion and relaxation. The changes in x are associated with big 
standard deviations due to small movements and tracking error. 



Table 1. Demographic and echocardiographic characteristics of 
subjects 



Parameters 


Mean 


±SD 


Body mass index (kg/rn) 


26.2 


± 2.0 


Heart rate (bpm) 


75 


± 5 


Diastolic pressure (mmHg) 


72 


± 14 


Systolic pressure (mmHg) 


111 


± 13 


Diastolic LV posterior wall (mm) 


11.9 


± 1.5 


Systolic LV posterior wall (mm) 


16.4 


±2.9 


Interventricularseptalend diastolic dimension 


12.0 


± 1.8 


(mm) 






Interventricularseptalend systolic dimension 


15.7 


± 3.0 


(mm) 






Fractional shortening (%) 


41.12 


±9.54 


Ejection fraction (%) 


62.92 


±9.82 


Diastolic LV large diameter (mm) 


79.5 


± 5.5 


Systolic LV large diameter (mm) 


64.0 


± 5.0 



LV: left ventricular 



Changes of z coordinates in both base and mid-levels are consis- 
tent with systolic and diastolic ventricular phases but are oppo- 
site to the expected positive direction of z in the apex level. 
Changes of 2 coordinates in the three base, mid and apex levels 
suggest that the septal wall contracts during the cardiac cycle, 
while moving in the direction of negative z (the direction of the 
apex). The average rates of temporal changes of X, 6 and O dur- 
ing three cardiac cycles for the mid-ventricular septal wall of a 
healthy individual are presented in Fig. 5 for three base (Fig. 
6 A), mid (Fig. 6B), and apex (Fig. 6C) levels in degrees. 

Considering the positivity contract of counter-clockwise direc- 
tion forms the apex point of view for 0, and negativity contract 
form perpendicular to the plane (_LZ, Z) point of view for X, the 
counter-clockwise rotation of the apex and clockwise rotation of 
the base are in agreement with the results of previous stud- 

3)6)7)11)15) non-definition of a unique direction for O, 

its magnitude is shown. 

The averages and standard deviations of Xm, 6m, and Om 
during three cardiac cycles for the interventricularseptal wall 
in healthy individuals in this study are presented in Fig. 7, for 
the three base, mid and apex levels in degrees. 

Considering Fig. 7, the bigger magnitude of Om in all three 
base, mid and apex levels compared to the other two angles is 
clearly obvious. 0m is smaller than the other two angles in 
both base and mid-levels which suggests a large septal wall 
movement in the longitudinal direction, compared to the short- 
axis direction. 0m was obtained bigger than the Xm angle in 
the apex level, and considering its large standard deviation, it 
requires fijrther investigation. 

After studying the normality of data (Shapiro- Wilk test), 
the correlations between Xm and Om, and also angle 0m in the 
three base, mid and apex levels were analyzed, using Pearson's 
analysis of correlation (Fig. 8). A highly significant correlation 
was observed between the two 0m and Om in the apex level (r 
= 0.82, p = 0.03). The results showed a significant correlation 
between Xm and Om in the base level of the interventricular- 
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septal wall (r = 0.57,p = 0.03) and in the mid-level of the in- 
terventricularseptal wall (r = 0.79,/^ = 0.01). 

The regression function between 0m and Om angles in the 
apex level was estimated as Om = 0.87 x 0m + 5.57. The regres- 
sion function between Xm and Om angles was estimated as Om 
= 0.35 X Xm + 11.09 in the base level of the interventricularsep- 
tal wall and as Om = 1.46 x Xm - 1.62 in the mid-level of the 
interventricularseptal wall. There was no significant correlation 
between 0m and Om angles, in both base and mid-levels. Om 
and Xm angles had no significant correlation in apex level. 

Discussion 

Quantification of the LV motion were primarily based on 
simplified geometries such as cylinder and truncated ellipsoid 
by considering the isotropic myofiber structure.^ ^"^^^ The geo- 
metric model of ellipsoid with isotropic myofiber structure was 
used in the current study, which is the closest simple geomet- 
ric model to the actual ventricular geometry. The LV's motion 
has been estimated by applying differential and area matching 




Fig. 7. Xm, 9m, and Om during three cardiac cycles for the interventri- 
cularseptal wall in healthy individuals for three base, mid and apex 
levels, in degrees. 



algorithms on sequential 2 -dimensional echocardiography or 
MRI images. The rotation angles for different segments 
were estimated using short-axis view images in two dimen- 
sions, and the mean of angles was considered as the rotation 
angle, for each level. The torsion angle (difference of angles in 
two base and apex levels),^ is a general parameter and de- 
pends on the changes in the whole ventricle. Since ventricle 
motion takes place in circumferential, radial and longitudinal 
directions, and because of the possibility of change in pattern 
of the ventricle motion inonly a particular direction, the rota- 
tion and torsion angles may not be sensitive enough in differ- 
ent physiological conditions. In order to have a parameter that 
includes all three directions, and reflects their changes, O was 
defined and estimated. A block matching algorithm was used. 
Factors such as changes in the shape of tissue, image noise, 
speckle non-correlation, off-the-screen motion, and artifacts 
caused by speckle motion, have restricted the use of block 
matching in ultrasonic images. 

Results of correlation between Om with 0m and Xm showed a 
significant correlation between the Om and 0m of apex level. 
Also Om and Xm were significantly correlated in both base and 
mid-levels. The non-correlation between the Om and 0m in 
both base and mid levels, and also the non-correlation be- 
tween Om and Xm in apex level could be due to specific orien- 
tation of fibers in those regions. According to Torrent-Gausp 
model, the fibers in base and mid regions of septum have lon- 
gitudinal orientation and the fibers in its apex segment have 
circumferential orientation.^^ Therefore, motion in base and 
mid levels of septal wall will be affected by longitudinal fibers 
and the X has a bigger contribution in O. The motion of sep- 
tal wall is also affected by circumferential fibers in apex level; 
therefore 0 has a bigger contribution in O. 

Gustafsson et al.^^^ reported the magnitude of 0 as -8.8 ± 
3.5° in the base, -2.3 ± 3.0° in the mid, and 11.1 ± 4.4° in the 
apex levels. 0m of the interventricular septum in the current 
study was estimated as -10.61 ± 3.38° in the base, -7.78 ± 
2.96° in the mid, and 10.37 + 5.48° in the apex segments. 
The range and sign of angles in both base and apex levels are 
in agreement with Gustafsson et al.^^^ A great difference is ob- 
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served in mid-level that isprobably associated with the lack of 
an accurate definition of mid-level. During the only study 
that was conducted in long-axis view, Arab-Baferani et al/^^ 
estimated Xm as 18.51 ± 3.41° in the base, and 14.74 ± 2.91° 
in the apex levels. Xm of the interventricular septum was esti- 
mated as 15.11 ± 3.30° in the base, 16.72 ± 2.66° in the 
mid, and 8.79 + 5.82° in the apex segments, respectively. It 
seems that the difference for apex level is associated with the 
lack of an accurate definition of this level, and tracking error, 
due to low quality of images. It should be noted that no stan- 
dard value has yet been defined for Xm and 0m. 

The magnitude of Om was estimated as 16.33 ± 3.01° in 
the base, 22.77 ± 4.95° in the mid, and 14.60 ± 5.81° in the 
apex levels. Urbano Moral et al.^^^ reported the average rota- 
tion in the base and apex levels as -2.6 ± 1.6° and +8.2 ± 3.6°, 
respectively using the speckle tracking of 3 -dimensional echo- 
cardiography images. This study was carried out using an Ar- 
tida 4D system (Toshiba Medical Systems, Tustin, CA, USA) 
that according to the instructions of its manufacturer, measures 
the rotation in two base and apex levels. Thus, the obtained re- 
sults cannot be compared with the results of this study. 

Due to the low temporal resolution of 3 -dimensional imag- 
es during echocardiography in the current study, and the im- 
possibility of saving raw 3 -dimensional images, the combina- 
tion of 2 -dimensional images was used. While having the 
advantage of high temporal resolution, this approach does how- 
ever increase the error for finding the location of the desired seg- 
ment in images of longitudinal view, and their corresponding 
segment in short-axis view, and suffers from imaging in two 
different times. Thus, using real-time 3 -dimensional images 
can be an effective way to remove some of these factors. It is 
suggested that such angle is measured using real-time 3 -di- 
mensional images and the capability of O to identify the 
physiological conditions of healthy and unhealthy hearts com- 
pared to the other two angles, will be examined. 

In conclusion, due to the motion of the wall of the LV in 
three dimensions, and the non-isotropic structure of myofi- 
bers, it is possible that in different physiological conditions, 
any measurable change only happens during the motion of 
myocardial muscle in one of the longitudinal, circumferential 
or radialdirictions. Therefore, a parameter dependent on the 
motion in three directions can be helpful. So, O for myocardi- 
al muscle was defined using the speckle tracking method of 
2 -dimensional echocardiography images in both short-axis 
and long-axis views and estimated in three base, mid and apex 
levels of the septal wall. The significant correlation between 
this angle and 0 and O show that this angle can be used to 
study the behavior of the myocardial muscle in different path- 
ological conditions. 
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